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Abstract 

Some physical implications involved in a new concept, termed the 
"quantum protectorate" (QP), are developed and discussed. This is 
done by considering the idea of quantum protectorate in the context 
of quantum theory of magnetism. It is suggested that the difficulties 
in the formulation of quantum theory of magnetism at the microscopic 
level, that are related to the choice of relevant models, can be un- 
derstood better in the light of the QP concept . We argue that the 
difficulties in the formulation of adequate microscopic models of elec- 
tron and magnetic properties of materials are intimately related to 
dual, itinerant and localized behaviour of electrons. We formulate a 
criterion of what basic picture describes best this dual behaviour. The 
main suggestion is that quasi-particle excitation spectra might provide 
distinctive signatures and good criteria for the appropriate choice of 
the relevant model. 
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1 Introduction 



It is well known that there are many branches of physics and chemistry 
where phenomena occur which cannot be described in the framework of 
interactions amongst a few particles |M. As a rule, these phenomena arise 
essentially from the cooperative behaviour of a large number of particles. 
Such many-body problems are of great interest not only because of the na- 
ture of phenomena themselves, but also because of the intrinsic difficulty 
in solving problems which involve interactions of many particles in terms 
of known Anderson statement that "more is different" Q. It is often dif- 
ficult to formulate a fully consistent and adequate microscopic theory of 
complex cooperative phenomena. In ref.Q, the authors invented an idea 
of a quantum protectorate, "a stable state of matter, whose generic low- 
energy properties are determined by a higher-organizing principle and noth- 
ing else" ||] . This idea brings into physics the concept that reminds the 
uncertainty relations of quantum mechanics . The notion of QP was in- 
troduced to unify some generic features of complex physical systems on 
different energy scales, and is a certain reformulation of the conservation 
laws and symmetry breaking concepts Q. As typical examples of QP, the 
crystalline state, the Landau fermi liquid, the state of matter represented 
by conventional metals and normal 3 He ( c/. || , JjJ ) , and the quantum Hall 
effect were considered. The sources of quantum protection in high-T c su- 
perconductivity and low-dimensional systems were discussed in refs.|8||-||To|l. 
According to Anderson ||, "the source of quantum protection is likely to 
be a collective state of the quantum field, in which the individual particles 
are sufficiently tightly coupled that elementary excitations no longer involve 
just a few particles, but are collective excitations of the whole system. As 
a result, macroscopic behaviour is mostly determined by overall conserva- 
tion laws" . In the same manner the concept of a spontaneous breakdown of 
symmetry enters through the observation that the symmetry of a physical 
system could be lower than the symmetry of the basic equations describing 
the system 0j,§|. This situation is encountered in non-relativistic statisti- 
cal mechanics. A typical example is provided by the formation of a crystal 
which is not invariant under all space translations, although the basic equa- 
tions of equilibrium mechanics are. In this article, I will attempt to relate 
the term of a quantum protectorate and the foundations of quantum theory 
of magnetism. I will not touch the low-dimensional systems that were dis- 
cussed already comprehensively in refs.||-|lC]. I concentrate on the problem 
of choosing the most adequate microscopic model of magnetism of materials 
and, in particular, related to the duality of localized and itinerant behaviour 
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of electrons where the microscopic theory meets the most serious difficul- 
ties. To justify this statement and to introduce all necessary notions that 
are relevant for the present discussion, we very briefly recall the basic facts 
of the microscopic approach to magnetism. 



2 Magnetic Degrees of Freedom 

The discussion in this paper is concentrated on the right definition of the 
fundamental "magnetic" degrees of freedom and their correct model de- 
scription for complex magnetic systems. We shall first describe the phe- 
nomenology of the magnetic materials to look at the physics involved. The 
problem of identification of the fundamental " magnetic" degrees of freedom 
in complex materials is rather nontrivial. Let us discuss briefly, to give 
a flavor only, the very intriguing problem of the electron dual behaviour. 
The existence and properties of localized and itinerant magnetism in insu- 
lators, metals, oxides and alloys and their interplay in complex materials is 
an interesting and not yet fully understood problem of quantum theory of 
magnetism) 11], [12]. The central problem of recent efforts is to investigate 



the interplay and competition of the insulating, metallic, superconducting, 
and heavy fermion behaviour versus the magnetic behaviour, especially in 
the vicinity of a transition to a magnetically ordered state. The behaviour 
and the true nature of the electronic and spin states and their quasi-particle 
dynamics are of central importance to the understanding of the physics of 
strongly correlated systems such as magnetism and metal-insulator transi- 
tion in metals and oxides, heavy fermion states , superconductivity and their 
competition with magnetism fj~3|| . The strongly correlated electron systems 
are systems in which electron correlations dominate. An important problem 
in understanding the physical behaviour of these systems was the connection 
between relevant underlying chemical, crystal and electronic structure, and 
the magnetic and transport properties which continue to be the subject of 
intensive debat es pjj j. Strongly correlated d and / electron systems are of 
special interest [fL5|] . In these materials electron correlation effects are es- 
sential and, moreover, their spectra are complex, i.e., have many branches. 
Importance of the studies on strongly correlated electron systems are con- 
cerned with a fundamental problem of electronic solid state theory, namely, 
with a tendency of 3(4)<i electrons in transition metals and compounds and 
4(5)/ electrons in rare-earth metals and compounds and alloys to exhibit 
both localized and delocalized behaviour |lH|,|]16[j. Many electronic and mag- 
netic features of these substances relate intimately to this dual behaviour 



2 



of the relevant electronic states. For example, there are some alloy systems 
in which radical changes in physical properties occur with relatively mod- 
est changes in chemical composition or structural perfection of the crystal 
lattice [15]. Due to competing interactions of comparable strength, more 
complex ground states than usually supposed may be realized. The strong 
correlation effects among electrons, which lead to the formation of the heavy 
fermion state take part to some extent in formation of a magnetically or- 
dered phase, and thus imply that the very delicate competition and interplay 
of interactions exist in these substances jlT]]. For most of the heavy fermion 
superconductors, cooperative magnetism, usually some kind of antiferro- 
magnetic ordering was observed in the "vicinity" of superconductivity. In 
the case of U-based compounds, the two phenomena, antiferromagnetism 
and superconductivity coexist on a microscopic scale, while they seem to 
compete with each other in the Ce-based systems) IS]. For a Kondo lat- 
tice system, the formation of a Neel state via the RKKY intersite interac- 
tion compete with the formation of a local Kondo singlet . Recent data 
for many heavy fermion Ce- or U-based compounds and alloys display a 
pronounced non-Fermi-liquid behaviour. A number of theoretical scenarios 
have been proposed and they can be broadly classified into two categories 
which deal with the localized and extended states of /-electrons. Of special 
interest is the unsolved controversial problem of the reduced magnetic mo- 
ment in Ce- and U-based alloys and the description of the heavy fermion 
state in the presence of the coexisting magnetic state. In other words, the 
main interest is in the understanding of the competition of intra-site (Kondo 
screening) and inter-site (RKKY exchange) interactions. Depending on the 
relative magnitudes of the Kondo and RKKY scales, materials with different 
characteristics are found which are classified as non-magnetic and magnetic 
concentrated Kondo systems. The latter, "Kondo magnets", are of main 
interest [15]. Furthermore, there are effects which have a very complicated 
and controversial origin. There are some experimental evidences that pe- 
culiar magnetism of some quasi-ternary heavy fermion alloys is not that of 
localized systems, but have some features of band magnetism. Thus, in 
addition to the pronounced non-Fermi-liquid effects in thermodynamic and 
transport properties, the outstanding problems include small magnetic mo- 
ments and possible transitions from a localized moment ordered phase to 
a kind of "heavy fermion band magnet" [IS] - pi]] . These features reflect 
the very delicate interplay and competition of interactions and changes in 
a chemical composition. As a rule, very little intuitive insight could be 
gained from this very complicated behaviour. The QP is an umbrella term 
for a theoretical approach which seems designed specifically to analyze such 
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problems. 

3 Microscopic Picture of Magnetism in Materials. 

In this Section we recall the foundations of the quantum theory of mag- 
netism in a sketchy form. Magnetism in materials such as iron and nickel 
results from the cooperative alignment of the microscopic magnetic moments 
of electrons in the material. The interactions between the microscopic mag- 
nets are described mathematically by the form of the Hamiltonian of the 
system. The Hamiltonian depends on some parameters, or coupling con- 
stants, which measure the strength of different kinds of interactions. The 
magnetization, which is measured experimentally, is related to the average 
or mean alignment of the microscopic magnets. It is clear that some of 
the parameters describing the transition to the magnetically ordered state 
do depend on the detailed nature of the forces between the microscopic 
magnetic moments. The strength of the interaction will be reflected in the 
critical temperature which is high if the aligning forces are strong and low 
if they are weak. In quantum theory of magnetism, the method of model 
Hamiltonians has proved to be very effective. Without exaggeration, one 
can say that the great advances in the physics of magnetic phenomena are 
to a considerable extent due to the use of very simplified and schematic 
model representations for the theoretical interpretation. 



3.1 Heisenberg Model 

The Heisenberg model is based on the assumption that the wave functions 
of magnetically active electrons in crystals differ little from the atomic or- 
bitals. The physical picture can be represented by a model in which the 
localized magnetic moments originating from ions with incomplete shells in- 
teract through a short-range interaction. Individual spin moments form a 
regular lattice. The model of a system of spins on a lattice is termed the 



Heisenberg ferromagnet[22] and establishes the origin of the coupling con- 
stant as the exchange energy. The Heisenberg ferromagnet in a magnetic 
field H is described by the Hamiltonian 

H = ~ E J (* " 3)$i$3 ~ WbH ]T Sf (1) 

ij i 

The coupling coefficient J(i — j) is the measure of the exchange interaction 
between spins at the lattice sites i and j and is defined usually to have 
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the property J(i - j = 0) =0. This constraint means that only the inter- 
exchange interactions are taken into account. The coupling, in principle, 
can be of a more general type (non-Heisenberg terms). For crystal lattices 
in which every ion is at the centre of symmetry, the exchange parameter has 
the property 

J(i-j) = J(j-i) 
We can rewrite then the Hamiltonian ((!]) as 

H=-J2J(i-j)(S?S* + S?Sr) (2) 

ij 

Here = S x ± iS y are the raising and lowering spin angular momentum 
operators. The complete set of spin commutation relations is 

[S+,Sr}_ = 2Sf% [S+,S-} + = 2S(S + l) - 2{Sff- 

[Sf, Sf ]_ = ±Sf8 ij - Sf = S{S + 1) - {Sf f - SrS+; 

{s +fs+i = Q) (sr) 2 S+i = o 

We omit the term of interaction of the spin with an external magnetic field 
for the brevity of notation. The statistical mechanical problem involving 
this Hamiltonian was not exactly solved, but many approximate solutions 
were obtained. 

To proceed further, it is important to note that for the isotropic Heisenberg 
model, the total z-component of spin Sf ot = J2i Sf is a constant of motion, 
i.e. 

[H,St ot ]=0 (3) 

There are cases when the total spin is not a constant of motion, as, for 
instance, for the Heisenberg model with the dipole terms added. 
Let us define the eigenstate \tpo > so that S^~\ipo >= for all lattice sites 
Ri. It is clear that > is a state in which all the spins are fully aligned 
and for which Sf |V>o >= S\ipo >. We also have 

J f = £e<**>J(i) = J_ 5 

i 

, where the reciprocal vectors k are defined by cyclic boundary conditions. 
Then we obtain 

H\^o >= " E J (* " i) s2 = -NS 2 J 
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Here N is the total number of ions in the crystal. So, for the isotropic 
Heisenberg ferromagnet, the ground state \ipo > has an energy —NS 2 Jq. 
The state l^o > corresponds to a total spin NS. 

Let us consider now the first excited state. This state can be constructed 
by creating one unit of spin deviation in the system. As a result, the total 
spin is NS — 1. The state 

is an eigenstate of H which corresponds to a single magnon of the energy 

E(q) = 2S(J - J q ) (4) 

Note that the role of translational symmetry, i.e. the regular lattice of spins, 
is essential, since the state > is constructed from the fully aligned state 
by decreasing the spin at each site and summing over all spins with the 
phase factor e lkR j (we consider the 3-dimensional case only). It is easy to 
verify that 

<^ k \Sl ot \^ k >=NS-l 



The above consideration was possible because we knew the exact ground 
state of the Hamiltonian . There are many models where this is not the 
case. For example, we do not know the exact ground state of a Heisenberg 
ferromagnet with dipolar forces and the ground state of the Heisenberg 
antifer r omagnet . 



3.2 Itinerant Electron Model 

E.Stoner has proposed an alternative, phenomenological band model of mag- 
netism of the transition metals in which the bands for electrons of differ- 
ent spins are shifted in energy in a way that is favourable to ferromag- 
netism. The band shift effect is a consequence of strong intra-atomic corre- 
lations. The itinerant-electron picture is the alternative conceptual picture 
for magnetism p3], pi]. It must be noted that the problem of antiferro- 



magnetism is a much more complicated subject [25]. The antiferromagnetic 
state is characterized by a spatially changing component of magnetization 
which varies in such a way that the net magnetization of the system is zero. 
The concept of antiferromagnetism of localized spins, which is based on the 
Heisenberg model and two-sublattice Neel ground state, is relatively well 
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founded contrary to the antiferromagnetism of delocalized or itinerant elec- 
trons . In relation to the duality of localized and itinerant electronic states, 
G.Wannier showed the importance of the description of the electronic states 
which reconcile the band and local (cell) concept as a matter of principle. 

3.3 Hubbard Model 

There are big difficulties in the description of the complicated problem of 
magnetism in a metal with the d band electrons which are really neither 
"local" nor "itinerant" in a full sense. The Hubbard model(l^] is in a certain 
sense an intermediate model (the narrow-band model) and takes into account 
the specific features of transition metals and their compounds by assuming 
that the d electrons form a band, but are subject to a strong Coulomb 
repulsion at one lattice site. The Hubbard Hamiltonian is of the form |2ql ,||27|l 

H = ^jtijAtT^a + f//2^n i(T n i _ a . (5) 

ijcr icr 

It includes the intra-atomic Coulomb repulsion U and the one-electron hop- 
ping energy tij. The electron correlation forces electrons to localize in the 
atomic orbitals which are modelled here by a complete and orthogonal set 
of the Wannier wave functions [<p(f — Rj)\- On the other hand, the kinetic 
energy is reduced when electrons are delocalized. The band energy of Bloch 
electrons is defined as follows: 

= N- 1 4 exp[ik{Ri - Rj] (6) 
k 

where N is the number of lattice sites. This conceptually simple model 
is mathematically very complicated p6|],||27||. The Pauli exclusion principle 
which does not allow two electrons of common spin to be at the same site, 
plays a crucial role. It can be shown, that under transformation RHR + , 
where R is the spin rotation operator 

i? = (g)exp(-^a i n) (7) 

3 

the Hubbard Hamiltonian is invariant under spin rotation, i.e., RHR + = H. 
Here <j> is the angle of rotation around the unitary axis n and a is the Pauli 
spin vector; symbol indicates a tensor product over all site subspaces. 
The summation over j extends to all sites. 
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The equivalent expression for the Hubbard model that manifests the prop- 
erty of rotational invariance explicitly can be obtained with the aid of the 
transformation 

oXXa^'^V ( 8 ) 
Z aa' 

Then the second term in (|5|) takes the following form 
As a result we get 

ijcr i 

The total z-component Sf ot commutes with Hubbard Hamiltonian and the 
relation (^) is valid. 

3.4 Multi-Band Models. Model with s — d Hybridization 

The Hubbard model is the single-band model. It is necessary, in principle, to 
take into account the multi-band structure, orbital degeneracy, interatomic 
effects and electron-phonon interaction. The band structure calculations 
and the experimental studies showed that for noble, transition and rare- 
earth metals the multi-band effects are essential. An important generaliza- 
tion of the single-band Hubbard model is the so-called model with s — d 



hybridization [28], [29]. For transition d metals, investigation of the energy 
band structure reveals that s — d hybridization processes play an impor- 
tant part. Thus, among the other generalizations of the Hubbard model 
that correspond more closely to the real situation in transition metals, the 
model with s — d hybridization serves as an important tool for analyzing of 
the multi-band effects. The system is described by a narrow d-like band, a 
broad s-like band and as-d mixing term coupling the two former terms. 
The model Hamiltonian reads 

H = H d + H s + H s ^ d (10) 

The Hamiltonian Hd of tight-binding electrons is the Hubbard model @. 

Hs = J24clc ka (11) 

ka 



s 



is the Hamiltonian of a broad s-like band of electrons. 

H s -d = V k{Aa a ko + a h°kv) ( 12 ) 
ku 

is the interaction term which represents a mixture of the d-band and s-band 
electrons. The model Hamiltonian ( |l0| ) can be interpreted also in terms of 
a series of Anderson impurities placed regularly in each site (the so-called 
periodic Anderson model ). The model (jlOj ) is rotationally invariant also. 



3.5 Spin-Fermion Model 

Many magnetic and electronic properties of rare-earth metals and com- 
pounds {e.g., magnetic semiconductors) can be interpreted in terms of a 
combined spin-fermion model |3(|,|3lJ that includes the interacting local- 
ized spin and itinerant charge subsystems. The concept of the s(d) — f 
model plays an important role in the quantum theory of magnetism, espe- 
cially the generalized d—f model, which describes the localized 4/(5/)-spins 
interacting with cZ-like tight-binding itinerant electrons and takes into con- 
sideration the electron-electron interaction. The total Hamiltonian of the 
model is given by 

H = H d + H d „ f (13) 

The Hamiltonian H d of tight-binding electrons is the Hubbard model (||). 
The term H d -f describes the interaction of the total 4/(5/)-spins with the 
spin density of the itinerant electrons 



H d-f = J °i&i = ~JN ^^E^-^L"^-" + Z <r S -q a ka a k+qo] 
i kq o 



where sign factor z a is given by 

Z a = (+,-) - CT = (T, |) 

and 



(14) 



5 — U 
-q 




a = + 
a = — 



In general the indirect exchange integral J strongly depends on the wave 
vectors J(k; k + q) having its maximum value at k = q = 0. We omit this 
dependence for the sake of brevity of notation. To describe the magnetic 
semiconductors the Heisenberg interaction term (|l]) should be added 33] 
( the resulting model is called the modified Zener model ). 
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These model Hamiltonians (and their simple modifications and combina- 
tions) are the most commonly used models in quantum theory of mag- 
netism. In our previous paper |l6|, where the detailed analysis of the neutron 
scattering experiments on magnetic transition metals and their alloys and 
compounds was made, it was concluded that at the level of low-energy hy- 
drodynamic excitations one cannot distinguish between the models. The 
reason for that is the spin-rotation symmetry. In terms of refs.||,[||], the 
spin waves ( collective waves of the order parameter ) are in a quantum 
protectorate precisely in this sense. I will argue below the latter statement 
more explicitly. 

4 Symmetry and Physics of Magnetism 

In many-body interacting systems, the symmetry is important in classifying 
different phases and understanding the phase transitions between them|4[,|||] 
. To implement the QP idea it is necessary to establish the symmetry prop- 
erties and corresponding conservation laws of the microscopic models of 
magnetism. The Goldstone theorem states that, in a system with broken 
continuous symmetry ( i.e., a system such that the ground state is not in- 
variant under the operations of a continuous unitary group whose generators 
commute with the Hamiltonian ), there exists a collective mode with fre- 
quency vanishing as the momentum goes to zero. For many-particle systems 
on a lattice, this statement needs a proper adaptation. In the above form, the 
Goldstone theorem is true only if the condensed and normal phases have the 
same translational properties. When translational symmetry is also broken, 
the Goldstone mode appears at zero frequency but at nonzero momentum, 
e.g., a crystal and a helical spin-density- wave (SDW) ordering. As has been 
noted, this present paper is an attempt to explain the physical implications 
involved in the concept of QP for quantum theory of magnetism. All the 
three models considered above, the Heisenberg, the Hubbard, and the spin- 
fermion model, are spin rotationally invariant, RHR + = H. The sponta- 
neous magnetization of the spin or fermion system on a lattice that possesses 
the spin rotational invariance, indicate on a broken symmetry effect, i.e., 
that the physical ground state is not an eigenstate of the time-independent 
generators of symmetry transformations on the original Hamiltonian of the 
system. As a consequence, there must exist an excitation mode, that is an 
analog of the Goldstone mode for the continuous case (referred to as " mass- 
less" particles). It was shown that both the models, the Heisenberg model 
and the band or itinerant electron model of a solid, are capable of describ- 
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ing the theory of spin waves for ferromagnetic insulators and metals [16]. 
In their paper [34], Herring and Kittel showed that in simple approxima- 
tions the spin waves can be described equally well in the framework of the 
model of localized spins or the model of itinerant electrons. Therefore the 
study of, for example, the temperature dependence of the average moment 
in magnetic transition metals in the framework of low-temperature spin- 
wave theory does not, as a rule, give any indications in favor of a particular 
model. Moreover, the itinerant electron model (as well as the localized spin 
model) is capable of accounting for the exchange stiffness determining the 
properties of the transition region, known as the Bloch wall, which separates 
adjacent ferromagnetic domains with different directions of magnetization. 
The spin-wave stiffness constant D is defined so that the energy of a spin 
wave with a small wave vector q is E ~ Dq 2 . To characterize the dynamic 
behaviour of the magnetic systems in terms of the quantum many-body the- 



ory, the generalized spin susceptibility (GSS) is a very useful tool|35|]. The 
GSS is defined by 

X(q,uj) = J dt« S~(t),St q » exp(-iwt) (15) 

For the Hubbard model = a J Off. This GSS satisfies the important sum 
rule 

Imx{q,oj)duj = 7r(nj, — 7ij) = — 2ir < S z > (16) 



It is possible to check thatfl6|] 

X (g» = _ 2<5 *> + 4{*(M - -([Q-,S+ ])} (17) 
u or q H y 

Here the following notation was used for qQ~ = [S~,H] and *ff(q,u) =« 
Qq\Q-q »w It is clear from (|l7|) that for q = the GSS ( |l5| ) contains only 
the first term corresponding to the spin-wave pole for q = which exhausts 
the sum rule (jl6|) . For small q, due to the continuation principle, the GSS 
x(q, w) must be dominated by the spin wave pole with the energy 

" = Dq 2 = —L—{q([Q-, <?+ ]) - q 2 lim liir (18) 

2 < b z > w^O ij^O 

This result is the direct consequence of the spin rotational invariance and is 
valid for all the three models considered above. 
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5 Spin Quasiparticle Dynamics 



In this Section, to make the discussion more concrete and to illustrate the 
nature of spin excitations in the above described models, let us consider 
the generalized spin susceptibility (GSS), which measures the response of 
"magnetic" degrees of freedom to an external perturbation [ 35 1. The GSS is 



expressed in terms of the double-time thermal GF of spin variables [22] [12], 
that is defined as 

x(g;t-t') =« S+(t),Sz q (t') »= -iO{t-t') < [S+(t),SZ q (t')U >= 

r+oo 

1/27T / duj exp(— iu)t)x(q, ui) (19) 



The poles of the GSS determine the energy spectra of the excitations in the 
system. The explicit expressions for the poles are strongly dependent on the 



model used for the system and the character of approximations [16], [35]. 
The next step in description of the spin quasiparticle dynamics is to write 
down the equation of motion for the GF. Our attention is focused on the spin 
dynamics of the models. To describe self-consistently the spin dynamics of 
the models one should take into account the full algebra of relevant operators 
of the suitable "spin modes", which are appropriate for the case. 

5.1 Spin Dynamics of the Hubbard Model 

Theoretical calculations of the GSS in transition 3d metals have been largely 
based on the single-band Hubbard Hamiltonian]^]. The GSS for this case 
reads 

X(g,w) =« (T+\aZ q » u (20) 

Here 

a k = a fcT afc +Pi ; a k =Y1 a H a k+ P 1 
p p 

The result of the RPA calculation 1 35] has the following form 

«(,,»)=«,>:,»„= I ^ (2i) 



where 



X %,u) = N-ij: Wfc J" nfc y (22) 

fe w + 4+ q - 4 - A 

A = Tf^2( n kl - "fct) ( 23 ) 
k 
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Table 1: EXPERIMENTAL DATA for TRANSITION METALS 



Element \ Data 


T 


D meVA 2 




A eV 


Qmax 


Fe 


1043 K 


280 


2.177 






Co 


1403 K 


510 


1.707 


0.91 




Ni 


631 K 


433 


0.583 


0.5 ±0.1 


0.8.4- 1 


MnSi 


30 K 


52 


0.4 







The excitation spectrum of the Hubbard model determined by the poles of 



susceptibility (22) is shown schematically in fig.l. The experimental data 
for three typical magnetic material are listed in Table 1. Note, that typically 
Qmax < 0-75k F . 

5.2 Spin Dynamics of the Spin-Fermion Model 

When the goal is to describe self-consistently the quasiparticle dynamics of 
two interacting subsystems the situation is more complicated. For the spin- 



fermion model (14) the relevant algebra of operators should be described 

by the 'spinor' (-.*) ("relevant degrees of freedom" ) |^lj . Once this has been 

done, one should introduce the generalized matrix spin susceptibility of the 

form _ 

f«S^\S_ k » «SZ\a_ k »\ . . . 

+ la- + - = x(k,oj) (24) 

The spectrum of quasiparticle excitations without damping follows from the 

poles of the generalized mean-field susceptibility. 

Let us write down explicitly the first matrix element xV 

<< SJ \S „ >> u - 



'' ' '' ' -o - JN- 1 ^ - n ; ) + 2J 2 N-V 2 < 5 2 > (1 - Uxf^Xo 

(25) 

X ^k^) = N-^ {np+ki ~ np]) (26) 



where 



U). 



p 



' P)fc = (w + ^-4 +fc -A) (27) 



A = 2JN- 1 / 2 <S^> -C/A^Vt - ni ) 
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This result can be considered as reasonable approximation for description of 
the dynamics of localized spins in heavy rare-earth metals like Gd. (c.f. 



The magnetic excitation spectrum that follows from the GF (|24p consists of 
three branches - the acoustic spin wave, the optic spin wave and the Stoner 



continuum plf . In the hydrodynamic limit, q — > 0, lo — > the GF (24) can 
be written as 

« S q \S_ q » - uj _ m (28) 
where the acoustic spin wave energies are given by 

1/2 £ fc (n fcT + n ki )(q^) 2 ei + (2A)" 1 £ fc (n fcT " n nWM^ 



E{q) ~ ~ 2iVV2 < SI > +(n T - n] ) 

(29) 

and 

<5g>=<5g>[l + 2 ^;"j» > ]-' (30) 

In GMF approximation the density of itinerant electrons ( and the band 
splitting A) can be evaluated by solving the equation 

= E < "far ** >= E[ ex P(^ e fe +C/iV ~ ln --- JAr ~ 1/2 < 5 o > -e^+ir 1 
fc fc 

(31) 

Hence, the stiffness constant D can be expressed by the parameters of the 
Hamiltonian (13). 

The spectrum of the Stoner excitations is given by [31] 

E s \q) = e d k+q -ef + A (32) 

If we consider the optical spin wave branch then by direct calculation one 
can easily show that 



E opt {q) = E° opt + D(UE opt /JA - l)q 2 



E» opt = J(ni-n i )+2J<S%> (33) 

From the equation (33) one also finds the GF of itinerant spin density in 
the generalized mean field approximation 

«<r+lcr- »°- jfffog) ( 34) 
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5.3 Spin Dynamics of the Multi-Band Model 

Now let us calculate the GSS for the Hamiltonian (|l0|). In general, one 
should introduce the generalized matrix spin susceptibility of the form 

V«S+|f7_ ? » «S q \sZ q »J 



Here 



q 



Let us consider for brevity the calculation of the Green function << a q |o"_„ >>. 
According to ref.[]35|], the object now is to calculate the Green function 
<< 9 k {q) = a\, +q ^a k ^\aZ q >> w . In the random phase approximation (RPA), 
the equations of motion for the relevant Green functions are reduced to the 
closed form 

(u + ef (k + q) - ef (k)) « 9 k (q)\a+ q » w = (n k+q[ - n k ^)A{q,u) (36) 

-V k+q « c{ +q[ a k] \at q >> w +V k « a\ +ql a k1 \aZ q >> w 
(u - ef(k) + e s k+q ) « c\ +ql a^\(T± q » w =< A+ q i a k+qi > a (q^) (37) 
+V k « cl +ql c k ^\ai q » w -V k+q « 9 k (q)\aZ q » w 
(w - ef (A: + q)- e%) « a\ +ql c k i;\ai q » w =< 4 T c fcT > A(q,u) (38) 
+V k « a\ +ql a k] \at q >> w -V k+q « c{ +ql c k ^\al q » UJ 

(w + e s k+q - e|) << c\ +qi c k] \at q >> w = (39) 
+V k « cl +ql a k] \at q »u> ~V k+q « a\ +ql c kr \at q >> w 

Here the following definitions were introduced 

4{k) = 4 + ¥rY,< a L a p*> ( 4 °) 

V 

A{q,u) = 1 - — « a q \at q » u 

To truncate the hierarchy of Green functions equations ( |36| ) - (^) the RPA 
linearization was used 

[9 k (q), H d ]_ ~ {e d k - e d k+q )9 k {q) + A9 k (q) (41) 

~^J2(< a k+ql a k+ql > ~ < 4t°*T >)9p(q) 
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[a\ +ql c kh H d ]^ ~ -e d k+q a\ +ql c k ^ - 



X(?, W ) =« a-\at q » u = - \ T J£!. A (42) 



Now, we will use these equations to determine the spin susceptibility of d- 
electron subsystem in the random phase approximation. It can be shown 
that 

X MF {qM 
l-U X MF (q,cj) 

We introduced here the notation x MF (q, UJ ) for the mean field susceptibility 

to distinguish it from the x^iq-, 1 ^) (|22|) . 

The expression for the X (?j w ) i s °f the form 

X MF {q,u) = ^£{(n fe+gi -n kr )[-\V k \ 2 ((u; + ef(k) + e s k+q ) (43) 
k 

+ (u, + ef(fc + g)-e|)) 
+(" + 4+ q ~ 4)(" + ef ( fe + " e k)( w " + e fc+?)] 

^ < 4+,i°*+9i > ( w + e f ( fc + 1) - 4)}}R~ 1 

where 

= {-|VfcI 2 ((o; + ef (fc + g) - ef (fc))(u; + ef (fe + g) - e|) (44) 
+( W _ ^(fc) + e | +9 )(a; + e% +q - 4) + (a> + ef (A + q) - e}(fc))( W - ef(fe) + e| +? ) 

+(w + ef(& + g) - e£)(w + e s k+q - e%j) 
+( W + e f(fc + g) - e}(fc))( W - ef(fe) + e| +9 )(u, + ef (fc + q) - e£)(u, + e£ +(? - e|)} 

Note, that if V k = then, x MF (<? 

, w) is reduced precisely to x O (<Z>k-0 (H)- 
The spectrum of quasiparticle excitations corresponds to the poles of the spin 
susceptibility (|22j) ; it corresponds to the spin- wave modes and to the Stoner- 
like spin-flip modes. Let us discuss first the question about the existence of 
a spin- wave pole among the set of poles of the susceptibility (p2|) . If we set 
q = in the secular equation for poles becomes 

1 = ^ E(Kl - n k1 )[-\V k \\2u - A) (45) 

k 

+u{uJ + e d {k) _ 4){uJ _ e d {k) + 4)] 
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-u[V k < 4 T c fcT > (w - ef (*;) + e|) + V* < 4 ; o fci > (w + ef (fc) - e|)]} 

(-|T4| 2 (2 W + A) 2 + 

w(w _ ejffc) + + e f (fc) - e |)(o; + 4 - e s k )) ' 

which is satisfied if u = 0. It follows from general considerations of Section 
4 that when the wave length of a spin wave is very long (hydrodynamic limit 
), its energy E(q) must be related to the wave number q by E{q) = Dq 2 . 
Thus the solution for the equation 

l = U X MF {q^) (46) 

exists which has the property linig^o E(q) = and this solution corresponds 
to a spin-wave excitation in the multiband model with s — d hybridization 



(f42p. Thus we derived a formula (42) for the dynamic spin susceptibility 
x(q,u)) in RPA and shown, that it can be calculated in terms of the mean 
field spin susceptibility X kQ^) by analogy with the single-band Hubbard 
model. 

Let us consider the poles of the x MF {l^ 00 )- It is instructive to remark that 
the Hamiltonian (|l^) can be rewritten in the mean field approximation as 

H MF = 4(kW ka a ka + J2 4cLc ka + J2 Vk{A a a ka + a^c fc(T ) (47) 

ka ka k<r 

The Hamiltonian ( f46| ) can be diagonalized by the Bogoliubov (u, ^-transformation 

a ka = UkoO-ko + v kcr (3 kcr ; c kcr = u ka j3 ka — v ku a ka (48) 
The result of diagonalization is 

H MF = Y,{uik*°L a ** + ^2kAjka) (49) 



ka 



where 



l/2[(ei(k) + e s k )±J(ei(k)-e° k )2 + 4\V k 



(51) 
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Then we find 



MF ( \ 1 V^r 2 2 ( n fct n k+qi) , eo> 

The present consideration shows that for the correlated model with s — d 
hybridization the spectrum of spin quasiparticle excitations is modified in 
comparison with the single-band Hubbard model. 



6 Quasiparticle Excitation Spectra and Neutron 
Scattering 

The investigation of the spectrum of magnetic excitations of transition and 
rare-earth metals and their compounds is of great interest for refining our 
theoretical model representations about the nature of magnetism. Experi- 
ments that probe the quasi-particle states could shed new light on the funda- 
mental aspects of the physics of magnetism. The most direct and convenient 
method of experimental study of the spectrum of magnetic excitations is the 
method of inelastic scattering of thermal neutrons . It is known experimen- 
tally that the spin wave scattering of slow neutrons in transition metals and 
compounds can be described on the basis of the Heisenberg model. On the 
other hand, the mean magnetic moments of the ions in solids differ appre- 
ciably from the atomic values and are often fractional. The main statement 
of the present consideration is that the excitation spectrum of the Hubbard 
model and some of its modifications is of considerable interest from the point 
of view of the choice of the relevant microscopic model. Let us consider the 
neutron scattering cross section which is proportional to the imaginary part 
of the GSS|5| 

d 2 a / 7e 2 ..o,— l.A;', 9, 

F(q)\ — )-(l + q 2 z ) 53 



dQduj ^m e c 2 ' 2 k 

[{N(u) + l)Imx(-q,w)+N(-u)Imx(q,w)] 

Here N(E(k)) is the Bose distribution function N(E(k)) = [exp(E(k)(3) - 
To calculate the cross section (53), we obtain from ( (42| ) the imaginary 
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part of the susceptibility, namely 



Imxiq u) = Im X MF ( q ,u;) 

The spin wave pole occurs where Imx MF (q, u>) tends to zero |35j| . In this 
case, we can in (|54l) take the limit Imx MF (q,u) — > so that 



UIm X (q, w) ~ -n5[l - URe X MF (q, u)] (55) 

but 

1 - URe X MF {q -» 0, cu -> 0) ~ b _1 (a; - £(g)) (56) 

and thus 

Jmx(g -> 0, u -> 0) ~ -tt^w - £(<?)) (57) 

Here 6 is a certain constant, which can be numerically calculated and E(q) 
is the acoustic spin wave pole E(q — > 0) = 0. 

Turning now to the calculation of the cross section (|53|), we obtain the 
following result 

^~(£)Vc«>Pc^a + «>4 « 58 » 

X)[^(£?(p))«y(w + S(p)) + (^(^(p)) + l)6(u - E{p))\ 
v 

According to formula (|58|), the cross section for the acoustic spin wave scat- 
tering will be identical for the Heisenberg and Hubbard (single-band and 
multiband) model. So, at the level of low-energy, hydrodynamic excita- 
tions one cannot distinguish between the models. However, for the Hubbard 
model, the poles of the GSS will contain, in addition to acoustic spin-wave 
pole, the continuum of the Stoner excitations E st {q) = e^+q — e g + A, as 
is shown on fig. |l[ The spectra of the spin-fermion model and multiorbital 
(multi-band) Hubbard model are shown for comparison. 
The cross section (p8|) does not include the contribution arising from the 
scattering by Stoner excitations, i.e. that determined by x MF (q^ UJ )- It was 



shown in paper [ 16 1 that in a single-band Hubbard model of transition metal 
in the limit when the wave vector of the elementary excitations goes to zero, 
the acoustic spin-wave mode dominates the inelastic neutron scattering, and 
the contribution to the cross section due to Stoner-mode scattering goes to 
zero. It was shown that the Stoner-mode scattering intensity does not be- 
come comparable to the spin- wave scattering intensity until q — 0-9<7ma:c 
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q max 




q 



Figure 1: Schematic form of the excitation spectra of the four microscopic 
models of magnetism: a) upper-left, the Heisenberg model; b) upper-right, 
the Hubbard model; c) down-left, the modified Zener (spin-fermion) model; 
d) down-right, the multiband Hubbard model. 

(see fig.l). Here q ma x is the value of q when the spin wave enters the contin- 
uum. For large values of q and u the energy gap A for spin flipping Stoner 
excitations may be overcome. In this case 

Imx(q,u) ~ Imx MF {q,uj) (59) 

From (^) we obtain for Imx MF {q,u) the result 

Im X MF (q, UJ ) = Jj- ^2{ u l+ql u l] ( n k^ ~ n% +gl )6(uj + u lk+ql - u^) (60) 

k 

+ V k+gl4^(n^ - n l+qi)^ + U2k+ql ~ W 2 fc T ) 

+ u k+ql v h( n k\ - n^ +ql )5(uj + ui k+ql - w 2 fct) + 

V k+ql U h( n k^ - ni +ql )5(uj + U 2 k+ql - OJlk])} 
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Now it follows from (^) that Imx MF (q,u) is non-zero only for values of 
the energies equal to the energies of the Stoner-type excitations 

(61) 



Ef\ q ) 


= u\k<[ 


- Ulk- 




El\q) 


= ^2fct 


~ U 2 k- 


Hi 


El\q) 


= ^2k\ 


- Wlfc. 


Hi 


Ef\q) 


= ^lfct 


- U 2 k- 


Hi 



With (J60j) and (|6_y) we obtain 



^ \ F (q)\ 2 (i)i-( 1 + ( iz)- (62) 



[(JV(w) + l)Im X MF (-q,u) + N{-ui)Im X MF {q,u>)] 

Although for the single-band model the Stoner-mode scattering cross sec- 
tion remains relatively small until q is fairly close to q m ax, it can be shown 



( see [16 1) that in the multiband models the Stoner-mode cross section may 
become reasonably large for much smaller scattering vector. 
The essential result of the present consideration is the calculation of the 
GSS for the model with s — d hybridization which is more realistic for tran- 
sition metals than the single-band Hubbard model. The present qualitative 
treatment shows that a two-band picture of inelastic neutron scattering is 
modified in comparison with the single-band Hubbard model. We have 
found that the long-wave-length acoustic spin-wave excitations should exist 
in this model and that in the limit (lim^^o nm <3-*o)) the acoustic spin- wave 
mode dominates the inelastic neutron scattering. The spin-wave part of the 
cross section is renormalized only quantitatively. The cross section due to 
Stoner-mode scattering is qualitatively modified because of occuring of the 
four intersecting Stoner-type sub-bands which may lead to the modifica- 
tion of the spin wave intensity fall off with increasing energy transfer. The 
intersection point q max can be essentially renormalized. 



7 Conclusions 

In summary, in this article, the logic of an approach to the quantum theory 
of magnetism based on the idea of the QP was described. There is an im- 
portant aspect of this consideration, which is seen to be the key principle 
for the interpretation of the spin quasiparticle dynamics of the microscopic 
models of magnetism. 

To summarize, the usefulness of the QP concept for physics of magnetism 
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derives from the following features. From our point of view , the clearest 
difference between the models is manifested in the spectrum of magnetic ex- 
citations. The model of correlated itinerant electrons and the spin-fermion 
model have more complicated spectra than the model of localized spins (see 
fig. ||) . Since the structure of the GSS and the form of its poles are de- 
termined by the choice of the model Hamiltonian of the system and the 
approximations made in its calculation, the results of neutron scattering 
experiments can be used to judge the adequacy of the microscopic models. 
However, it should be emphasized that to judge reliably the applicability of 
a particular model, it is necessary to measure the susceptibility (the cross 
section) at all points of the reciprocal space and for a wide interval of tem- 
peratures, which is not always permitted by the existing experimental tech- 
niques. Thus, further development of experimental facilities will provide a 
base for further refining of the theoretical models and conceptions about the 
nature of magnetism. In terms of ref.||], to judge which of the models is 
more suitable, it is necessary to escape the QP. This can be done by mea- 
surements in the high (q, u>) region, where (q 

^ Qmax i 

E~ A) . 

The following statements can now be made as to our analysis and its 
results. In this paper, we shown that quasiparticle dynamics of magnetic 
materials can be reasonably understood by using the simplified, but work- 
able models of interacting spins and electrons on a lattice in the light of 
the QP concept. The spectrum of magnetic excitations of the Hubbard 
model reflects the dual behaviour of the magnetically active electrons in 
transition metals and their compounds. The general properties of rota- 
tional invariance of the model Hamiltonians show that the presence of a 
spin-wave acoustic pole in the generalized magnetic susceptibility is a direct 
consequence of the rotational symmetry of the system. Thus, the acous- 
tic spin-wave branch reflects a certain degree of localization of the relevant 
electrons; the characteristic quantity D, which determines the spin wave 
stiffness, can be measured directly in neutron experiments. In contrast, in 
the simplified Stoner model of band ferromagnetism the acoustic spin-waves 
do not exist. There is a continuum of single-particle Stoner excitations only. 
The presence of the Stoner continuum for the spectrum of excitations of 
the Hubbard model is a manifestation of the derealization of the magnetic 
electrons. Since the Stoner excitations do not arise in the Heisenberg model, 
their direct detection and detailed investigation by means of neutron scat- 
tering is one of the most intriguing problems of the fundamental physics of 
magnetic state. Concerning the QP notion studied in the present paper, an 
important conclusion is that the inelastic neutron scattering experiments on 
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metallic magnets permit one to make the process of escaping the QP very 
descriptive. In this consideration, our main emphasis was put on the aspects 
important from the point of view of quantum theory of magnetism, namely, 
on the dual character of fundamental " magnetic degrees of freedom" . Gen- 
erally speaking, the fortunate circumstance in this discussion is the fact that 
besides the very general idea of QP also concrete practical tools are available 
in the physics of magnetism, and the combination of these two approaches 
is possible in the neutron scattering experiments ( for details see ref.[15]). 
The approach is very versatile since it uses the symmetry properties in the 
most ingenious fashion. By this consideration an attempt is made to link 
phenomenological and quantum theory of magnetism together more firmly, 
thus giving a better understanding of the latter. Finally, to clarify the con- 
cept of QP, we comment on somewhat resembling mathematical structures 
which are encountered when one tries to implement classical dynamic sym- 
metries in quantum field theory within these schemes one is trying to 
fit a classical description of particles endowed with internal structures, like 
spin. However, these analogies, as well as the elaboration of an adequate 
mathematical formalism for expression of the concept of QP need further 
studies. Further work is also necessary for the development of compact 
criteria appropriate for the QP occurrence in all applications. 
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